Classical to Quantum Transition of a Driven Nonlinear Nanomechanical Resonator 
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We seek the first indications that a nanoelectromechanical system (NEMS) is entering the quan- 
tum domain as its mass and temperature are decreased. We find them by studying the transition 
from classical to quantum behavior of a driven nonlinear Duffing resonator. Numerical solutions of 
the equations of motion, operating in the bistable regime of the resonator, demonstrate that the 
quantum Wigner function gradually deviates from the corresponding classical phase-space proba- 
bility density. These clear diff"erences that develop due to nonlinearity can serve as experimental 
evidence, in the near future, that NEMS resonators are entering the quantum domain. 
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The race to observe quantum-mechanical behavior in 
human-made nanoelectromechanical systems (NEMS) is 
bringing us closer than ever to testing the basic prin- 
ciples of quantum mechanics [l|, Q. With recent ex- 
periments coming within just an order of magnitude 
from the ability to observe quantum zero-point mo- 
tion [slpyjpl, ideas about the quantum-to-classical tran- 
sition [y, lZI, IsI may soon become experimentally accessi- 
ble, more than 70 years after Schrodinger described his 
famous cat paradox [9]. As nanomechanical resonators 
become smaller, their masses decrease and natural fre- 
quencies Q increase — exceeding IGHz in recent experi- 
ments [lol, El. For such frequencies it is sufficient to 
cool down to temperatures on the order of 50mK for the 
quantum energy HQ to be comparable to the thermal en- 
ergy ksT. Cooling down to such temperatures should 
allow one to observe truly quantum mechanical phenom- 
ena, such as resonances, oscillator number states, super- 
positions, and entanglement [l2|, [l^, [HI, , at least 
for macroscopic objects that are sufficiently isolated from 
their environment. 

In this Letter we seek the ffist experimental signa- 
tures, indicating that a mechanical object is entering the 
quantum domain as its mass and temperature are de- 
creased. Are there any clear quantum-mechanical cor- 
rections to classical dynamics that we can hope to ob- 
serve even before the fully-quantum regime is accessi- 
ble experimentally? To answer this question we perform 
two separate calculations on the same physical object — a 
driven nanomechanical resonator — viewing it once as a 
quantum-mechanical system and once as a classical sys- 
tem. With everything else being equal, this allows us to 
contrast the dynamics of a classical resonator with that 
of its quantum clone. We can then search for a regime 
in which quantum dynamics just begins to deviate away 
from classical dynamics, providing us with the experi- 
mental evidence we are looking for. 

We calculate the dynamics numerically. We start 
each quantum calculation with a coherent state, which 
is a minimal wave-packet centered about some point in 



phase space. We start the corresponding classical calcu- 
lation with an ensemble of initial conditions — typically 
of N = 10^ points — drawn from a Gaussian distribution 
in phase space that is identical to the initial quantum- 
mechanical probability density. We display the calcu- 
lated quantum dynamics in phase space using the quan- 
tum Wigner function 
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dx'e-^P''\x^x'\p{t)\x-x'), (1) 



where p{t) is the density operator, and compare it with 
the time-evolution of the classical phase-space density. 
We remind the reader that the Wigner function is not 
a true probability distribution as it may possess nega- 
tive values, particularly when the quantum state has no 
classical analog. Nevertheless, its governing equation of 
motion reduces to the classical Liouville equation upon 
formally setting h = 0, or whenever the potential is no 
more than quadratic. Moreover, it reduces to the quan- 
tum probability P(x, t) of observing the system at posi- 
tion X at time t upon integration over p, and vice versa. 

We perform our calculations for three qualitatively dif- 
ferent situations: (1) An isolated resonator with no cou- 
pling to the environment; (2) A resonator coupled to a 
heat bath at temperature T^nv = 0; and (3) A resonator 
coupled to a heat bath at a finite temperature T^nv > 0- 
In all cases we expect to find a regime in which the evolu- 
tion of quantum observables agrees with the correspond- 
ing classical averages at least up to the so-called "Ehren- 
fest time" pjl, [l8|, although it is generally difficult to 
strictly define this time scale [HI, [20]. For the isolated 
resonator, we do not expect to see a convergence of the 
two phas e-sp ace distributions, as the limit of ^ is 
singular 21, [22I]. When we couple the resonator to an 
environment, as in a real experiment, we do expect to 
see a gradual classical to quantum transition |22|, [2^, [ij] . 
We wish to study the details of this transition here. 

As mentioned above, the quantum dynamics of a co- 
herent state in a harmonic potential is essentially classi- 
cal. We therefore consider a driven Duffing resonator — 



2 




(a) Initial coherent state, t = 



(b) t = 30T 




(c) t = 156T 



(d) t = 300T 



FIG. 1: Isolated driven Duffing resonator with h — 0.2. 
Wigner functions (left) and classical distributions (right) of 
the initial minimal Gaussian wave packet and its evolution at 
three later times. A square of area h is shown at the bottom 
right corner of (a). 



a nonlinear resonator commonly observed in experi- 
25 , 26|, 27, 2^, whose Hamiltonian is given by 



ments 



H, 



sys 



+ - xF cos ut + (2) 



where the mass m of the resonator, and its natural fre- 
quency Q have been scaled to 1 by an appropriate choice 
of units. Thus, we measure the degree in which we 
approach the quantum domain by an increasing effec- 
tive value of as measured in the scaled units. As we 
shall see, quantum corrections clearly appear for values of 

:^ 0.1. We shall re-interpret these values later in terms 
of real experimental masses and frequencies. The remain- 
ing parameters in the Hamiltonian are the driving am- 
plitude F, the driving frequency cj, and the non-linearity 
parameter e. For the calculations presented here we take 
s = 0.01, F = 0.015 - 0.06, and cj = 1.016 - 1.02. For 
this choice of parameters the classical Duffing resonator 
is in the bist ability regime, where in the presence of dissi- 
pation it can oscillate at one of two different amplitudes, 
depending on its initial conditions. 

For the isolated quantum resonator we integrate the 
Schrodinger equation numerically by expanding the wave 
function in a truncated harmonic-oscillator energy basis 
|n). From the wave function we calculate the Wigner 
function ([T]), and compare it to the integrated classical 
trajectories of the corresponding Gaussian distribution of 
initial conditions. Figure [1] shows the two distributions, 
calculated for = 0.2, at t = and at three later times, 
measured in units of driving period T = 27r/a;. The 
Wigner functions are scaled to a power of 1 /4 to enhance 
weak features, where blue denotes positive values, and 
red denotes negative values. A square of area h is given 
at the bottom right corner of Fig. [TJ^a). 

We clearly see similarities between the quantum and 
the classical distributions, namely the positive backbone 



of the Wigner function which resembles the classical dis- 
tribution for short times, and the maxima of the two 
distributions, which are roughly in the same positions in 
phase space. Nevertheless, two differences prevent the 
distributions from approaching each other even in the 
limit of h ^ 0. These are the strong interference pat- 
tern in the Wigner function and the infinitely fine struc- 
ture that develops in the classical distribution. As noted 
by Berry [2l|, the only way to obtain a smooth tran- 
sition from classical to quantum dynamics in this case 
is to perform some averaging over a finite phase-space 
area AxAp taking into account the limited preci- 

sion of typical measurement devices. Such averaging will 
smooth out the delicate structure in the classical phase- 
space density and will cancel out the interference fringes 
in the quantum Wigner function, bringing the two into 
coincidence. 

The need for averaging is avoided once we consider 
the influence of an environment. For the classical res- 
onator we use a standard Langevin approach, adding a 
velocity dependent dissipative term — jx and a time de- 
pendent random force SF{t). The latter is assumed to 
be a ^-correlated Gaussian white noise, related to the 
dissipative term according to the fluctuation-dissipation 
theorem, thus defining a temperature Tenv The quan- 
tum resonator is coupled linearly to a bath of simple har- 
monic oscillators in thermal equilibrium at temperature 
Tenv This adds two terms to the quantum Hamiltonian 
([2j) — a Hamiltonian for the bath Hi^athi and a standard 
interaction Hamiltonian [15], taken to be of the Caldeira- 
Leggett [29] type in the rotating-wave approximation. 



Hint = ^{K.ibia^ + 



(3) 



where the are bilinear coupling constants, the bi are 
annihilation operators acting on the bath oscillators, and 
a is the annihilation operator of the Duffing resonator. 

We consider the reduced density operator psys of the 
Duffing resonator by tracing over the bath degrees of free- 
dom. By assuming the interaction to be weak, and by 
employing the Markov approximation which assumes the 
bath has no memory, we obtain a standard master equa- 



tion 



psys 



— [Hsys.Psys] ~ 2 ^ n){a^apsys + PsysCt^CL 



-2apsysa^) - -n{aa^f 



+ psysdCL^ - 2o) PsysO), (4) 



where ft = (g^^/^sTer^^ _ i^-i ^j^^ Bose-Einstein mean 
occupation, 7 = 27rg{rt)\tz{Q)\'^ is interpreted as the 
damping constant, and g{ft) and tv{ft) are the density 
of states of the bath and the coupling strength, respec- 
tively, both evaluated at the natural frequency Q = 1 of 
the resonator, because the nonlinearity is small. Instead 
of integrating (j3|) directly, we use the Monte- Carlo wave 
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(a) Initial coherent state, t = 



(b) t = lOT 



(c) t = SOT 



(d) t = 300T 




(a) Initial coherent state, t = 



(b) t = SOT 




(c) t = i8or 



(d) t = SOOT 



FIG. 2: As in Fig. [T] but for a driven Duffing resonator with 
h = 0.1, coupled to a heat bath at Tenv = 0, with damping 
constant 7 = 0.01. The stable solution towards which all 
initial classical points flow is encircled in (d). 



FIG. 3: As in Fig. [T] but for a driven Duffing resonator with 
h = 0.1, coupled to a heat bath at ksTenv = with 
damping constant 7 = 0.001. The initial state is centered on 
the separatrix so that classically it splits towards both of the 
stable solutions. 
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which is more efficient compu- 



function method 
tationahy. 

Figure [2] shows the calculated results for a Duffing res- 
onator with h = 0.1, coupled to a heat bath at tem- 
perature Tenv = 0- We take the initial state to be cen- 
tered around a point in phase space that under classi- 
cal dissipative dynamics flows towards the state of large 
amplitude oscillations. At short times we see a general 
positive outline of the Wigner function which is similar 
to the classical density, but it quickly deviates from the 
classical distribution. The uncertainly principle prevents 
it from shrinking to a point as in the classical dynamics. 
More importantly, we clearly observe that the Wigner 
distribution has substantial weight around the state of 
small-amplitude oscillations, which is inaccessible classi- 
cally for the chosen initial conditions at T^nv = 0. The 
quantum resonator can switch between the two stable 
states via tunneling or quantum activation [34, al- 
though at this point it is impossible for us to distinguish 
between these two processes. 

Figure [3] shows the calculated results for a Duffing res- 
onator with h = 0.1, coupled to a heat bath at a finite 
temperature ksTenv = 2hQ. This temperature is ob- 
tained by adding a fluctuating force while reducing the 
dissipation, thus requiring a longer time for the resonator 
to reach its final steady state. Here we choose the initial 
state to straddle the separatrix between regions in phase 
space that flow to the two stable states. The interference 
pattern that develops for short times within the general 
positive outline in the Wigner function is soon destroyed 
by decoherence, and becomes erratic in space and time. 
At long times, both distributions peak around the two 
stable states, nevertheless they differ significantly. The 
classical density is tightly localized around the two solu- 
tions with no overlap, indicating that T^nv is too small 
to induce thermal switching between the states, as was 



FIG. 4: Quantum (left) and classical (right) probability den- 
sities P{x) for the phase-space distributions in Fig.[3jd). The 
y axes are scaled differently for better visualization. 



observed in a recent experiment [28]. The Wigner func- 
tion, on the other hand, is spread out in phase space, 
indicating that h is sufficiently large for the quantum 
resonator to switch between the two states via tunnel- 
ing or quantum activation [3^5 IIH . This is demonstrated 
more clearly in Fig.[4l which shows the probability distri- 
butions P{x) at the steady state. We find that only for 
temperatures as high as ksT = 17 HQ does the classical 
phase-space distribution become as wide as the Wigner 
function is at ksT = 2hrt, as shown in Fig. [5l Thus, in 
a real experiment, evidence for quantum- mechanical dy- 
namics can be demonstrated as long as temperature and 
other sources of noise can be controlled to better than an 
order of magnitude. 

We have found the evidence we were seeking for a me- 
chanical system entering the quantum domain with ef- 
fective h = 0.1 and ksTenv = 2hQ. It stems from the 
ability of a quantum Duffing resonator to switch between 
the two stable states, having finite probability of being 
in between the two states, while the classical resonator 
cannot. We intend to study individual quantum trajec- 
tories of the resonator to determine whether switching 
takes place via tunneling or quantum activation. 

An effective :^ 0.1 implies that the area in phase 
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(a) Phase space densities, t = SOOT 





(b) Probability densities P(x) at t = SOOT 



FIG. 5: As in Figs. [3jd) andlH only that the classical calcu- 
lation (right) is for ksTenv — IT/iQ, yielding similar phase- 
space distributions (a) and probability densities (b). Both 
distributions in (a) are left unsealed. 



space where the dynamics takes place, is on the or- 
der of lOfi. This area is roughly XmaxPmax — mVtac^^ 
where ac — the critical oscillation amplitude for the on- 
set of bist ability — is proportional to d/^/Q [26], with d 
being the diameter of the resonator , and Q its quality fac- 
tor. In current NEMS resonators [ii, Q m 10"^^kg, 
^ lO^Hz, and ac — 10~^m, yielding an effective h 
of 10~^, which is well within the classical domain. Yet, 
if we consider suspended nanotubes [3^, Is^] , then opti- 
mistic values may give m :^ 10~^^kg, lO^Hz and 
ac — 10~^^m, for which the effective :^ 10, which is 
much better than needed, while requiring a temperature 
of 

Tenv — lOmK. We therefore believe that it should be 
possible to realize our calculations in a real experiment 
in the near future. 
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